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Abstract 

The higher order moments of the fluctuations for the thermodynamical systems in the presence 
of fields are investigated in the framework of a theoretical method. The method uses a generalized 
statistical ensemble consonant with the adequate expression for the internal energy. The appli- 
cations refer to the case of a system in magnetoquasistatic field. In the case of linear magnetic 
media one finds that for the description of the magnetic induction fluctuations the Gaussian ap- 
proximation is good enough. For nonlinear media the corresponding fluctuations are non-Gaussian, 
they having a non-null asymmetry. Additionally the respective fluctuations have characteristics of 
leptokurtic, mesokurtic and platykurtic type, depending on the value of the magnetic field strength 
comparatively with a scaling factor of the magnetization curve. 
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I. INTRODUCTION 



In our previous work ||T| we have presented a phenomenological approach of the fluctua- 
tions for the generahzed systems (i.e. thermodynamical systems in the presence of fields). In 
the respective work the fluctuations were evaluated only by second order numerical character- 
istics (correlations and moments). From a more general probabilistic perspective 0, ^, ^, ^, ^ 
the fluctuations in physical systems must be evaluated also by means of higher order numer- 
ical characteristics (higher order correlations and moments). In the present paper we aim 
to present an approach for the evaluation of the higher order correlations for the alluded 
generalized systems. 

For our aim in Sec. II we will present the general considerations and relations regarding 



the approached problems. Next one in Sec.|T| we proceed to evaluations of higher order 
moments for linear magnetic systems (in the presence of magnet oquasist at ic field), as well 
as for nonlinear magnetic media. The mentioned evaluations are finalized in some interesting 
relations able for comparations with adequate experimental results. 



II. THEORETICAL CONSIDERATIONS 

A. Statistical ensembles for generalized systems 

As in previous work (and above reminded) we consider a generalized system described 
by the set of extensive parameters {U, Xi, X2, ■ . . , X„, ¥±,¥2, . . . , ¥m), where U denotes the 
generalized internal energy, Xi{i = 1,2, .. . ,n) signify the usual extensive parameters (for 
systems in absence of fields), while ¥j {j = 1,2,... ,m) mean the extensive parameters 
associated with the fields. 

In the framework of fluctuation theory the thermodynamical quantities represent the 
mean (or expected) values of random variables. In the following we will note with A the 
mean value of the random quantity A. In such context we have U = E, where E denotes 
the energy regarded as random variable. 

We take the investigated system as a small part of a big, isolated ensemble. Generally 
speaking the probability density for the fluctuations of the quantities E, {Xj}"^-^ and {¥j}"^^^ 
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are of the form [0, ||, |[ 



w = exp ( -PE - J2 - Yl ^J-^i ) 

2=1 j = l 



(2) 



where Z is given by the normahzation condition of the probabihty. 

■■■ / r[dX,exp(-«,X,)r[dF,exp(-7,r,) / e'^^^dl 
i=i j=i -^r 

dr represents the elementary volume in phase space. 

The statistical integral Z is a function of the quantities /?, {ai}^^^ and {7j}^]^. For the 
identification of the physical signification of the parameters (3, ai and 7^ we note 



m 



i=i j=i 



n m 



(3dE + J2 + J2 (3) 

i=i j=i 



or 



m 



i=i j=i 

n m 

= /?df/ + ^«idX;+^7,dl^ (4) 

i=i j=i 

where we take into account that in fact U = E. 

From the thermodynamics of the generalized systems it is known [Tl|] that 

n m 

df/ = Td5 + ^£dX^+^^jdT^ (5) 

i=i j=i 

where S, T and denote respectively the entropy, temperature and the field dependent 
intensive parameters, while ipj signify the supplementary parameters due to the presence of 
the fields. 

From @ and @ one obtains: 

d I In Z + f3U + apQ + ^ 7^1^ 

\ i=l 3=1 

n m 

= l3TdS + 5^ («2 + 13^ dX, + Y, (7, + f^i^j) ^Yj (6) 
i=i 3=1 
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One observes that in this relation the left hand term is an exact differential. This fact impose 
the condition that the right hand term of the equation to be also an exact differential. Prom 
the mentioned condition it directly results that we must have the following relations: 

oci = -Pi = (8) 

7,- = -/3V',- = (9) 

In the above relations k denotes the Boltzmann's constant. 

We observe that the quantities {oLi}^^i are functions of the field dependent intensive 
parameters. 



B. The evaluation of the higher order correlations and moments 



The mean values of the fluctuating quantities can be evaluated through the relations: 



a(lnZ) 
d[5 



i = 1,2, . . . ,n 
j = 1,2,... ,m 



(10) 



d{\nZ) 



l3,Oll,l3 



(11) 



d{\nZ) 



1^3 



(12) 



The expressions of the type Ilil^"^*)'^' Iljl^^)'^^ called higher order correlations. By 
using the statistical sum Z as above introduced for some of the respective correlations one 
obtains the expressions: 



daadab da a 



a,b = 1,2,... , n 



(13) 



d^ad% Ola Oljja 



(14) 
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a = 1, 2, . . . , 



d^ilnZ) 

daadjb daa '"^ ' 6 = 12 



(15) 



6XaSXbSXc 



daadatdac doadab 



a,b,c = 1, 2, ... n (16) 



d^ad^bd^c d^ad^b 



d^jjadljjb 



a,b,c = 1, 2, ... m (17) 



a = 1, 2, ... n 



daad'jbd'jc daad'jb 



(18) 



<9^a<9V'6 6,c = 1,2, . . . ,m 

The formulae for the correlations of orders higher than 3 are generally more complicated. 
But the higher order moments can be obtained by means of the following recurrence formulae 



0: 



daa daa 



(19) 



{^YaY 



+1 



d 



dY 

i6Y,r - ni6Yar~^ 



As examples we give here the expressions of the moments {6Xa)^ and {6Ya)^. 



(SXa) 



dal 
d 



+ 3 



dXa 
daa 



kT^ I Xa + ?>{kT 



. dXg 

dia 



(20) 



(21) 



d^'Ya 



+ 3 



dll \d-fa 



dYa 



kT- 



d 



Ya + 3{kT 



Ma 



(22) 



d'ipaj \ dlpa. 

The fourth order moments are of interest for the evaluation of the so called excess coef- 
ficient: 



{SXa 



(23) 



(sxar 

which is an indicator of the deviation from the gaussian distribution 



5 



III. HIGHER ORDER MOMENTS FOR THE SYSTEMS SITUATED IN A 
MAGNETOQUASISTATIC FIELD 



A. Linear magnetic media 

Let us consider a uniformly magnetized continuous media, situated in a magnetoqua- 
sistatic field. The system is characterized by the extensive parameters {U,V, N,'B), where 
V, N and B denote respectively the volume, particle number and magnetic induction. In 



the case of linear magnetic systems the differential of the internal energy is given by [[I^, |Tl 



dU = TdS - pdV + CdN + VU dB (24) 

For the sake of brevity we omitted the mean symbol from above the parameters V, iV, B, 
but we will take into account that in fact in relation ( p^ and in the expressions of the 
moments appear the mean values of the respective quantities. 

The relations (0) and (|T6|) show that the moments associated to the usual thermody- 
namic quantities, i.e. volume V or particle number N, are functions of the parameters ^j, 
which depend on the fields constraints. 

For example, in the case B = const, for the volume V one obtains: 

dV 



{SVy = ~kT{ — ] (25) 



T,C,B 



= k'T^ { ^ ) (26) 



where 10 



piB = const.) = pB,N = p - in ■ B - Ih^P^ (27) 



N 



H signify the magnetic field strength, fi is the magnetic permeability and p — y. 

By using the properties of the Jacobians the relation (pSj) can be transformed as follows: 

diV, I T,B) , ^ diV, I T, B) diV, N, T, B) 



(svy = -kT ^jr ' ' = -kT 



C,r,B) 5(^,iV,T,B) 9(p,C,T,B) 



dN J 



-kT ^ — '-^^ 5 (28) 



J T,N,B J Ty,B 



In the above relation we take into account the condition — ( Mf] = i ^] ■ Here is 



the place to point out that the result ( pSj) is the same with the one obtained within the 
Gaussian approximation. 

Now let us evaluate the second and third order parameters of fluctuations for the magnetic 
induction B. For simplicity we consider that volume and particle number are constant. By 
using the relations (p^lf), (|17D and (|22|) we find: 



kTV fd'^B\ fkrV dfi 



The result (p9D is identical with the one obtained within the Gaussian approximation. 



We remark that in the case of linear magnetic media {6B)^ is null, because /i is independent 
on H. Additionally in this case the excess coefficient (|23|) is also null. These facts show that 
in the alluded case the Gaussian approximation is sufficient for a quantitative description of 
the fluctuations for B. 



B. Nonlinear magnetic media 

In the case of nonlinear magnetic media /i depends on H. Therefore the evaluation of 
the moments of orders higher than 2 become necessary. 

We approach such a case under the constraints when V = const, and TV = const. Then 
for the internal energy U we have: 

dU = TdS + 1/H ■ dB 

For the moments of 2, 3 and 4 order of B we obtain: 
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'VJ VdH^J 



+ 3 



Ty,N 



{6By 



(33) 



For finding the explicit expressions of {6B)'^, (55)^ and {SB)* it is necessary to know 
tlie expression of the function B = B{H). The most known such function is given by the 
Langevin equation: 



where 



B = ^oMs ( coth a ) + fioH 



kT 



(34) 



(35) 



Ms represents the saturation magnetization, /io is the vacuum permeabihty and m signify 
the magnetic moment of an individual molecule. 

By means of some simple mathematical operations one finds: 



{5By 



kTfiQ 
V 



kT 



1 



1 



sinh a 



1 



{5Bf 



2filrrrMs ( cosh a 



V sinh a 



(36) 



(37) 



{SBY 



2^^171 Ms j 3 sinh a — 3 cosh a 



sinh a 



+ 



fiQmMs f 1 



1 



kT \a? sinh^a. 
One of the most used dependence of B from H is given [|12| by the function 



B = fioMs 



1 — exp 



2a2 



(38) 



(39) 



where a is a scaling factor. 

The mentioned function imply a maximum for the differential permeability fid = dB / dH 
at the value H = a. 

By using the general formulas for the 2nd, 3rd and 4th order moments of the random 
variable B one obtains: 



-—2 fiokT r MsH 
(SB) =^^<j_^exp 



+ 1 



(40) 
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= ir ^^^^-3 exp 



In the end we wish to note the following observations: 



1. {5By change its sign in the point H = a, where the differential permeability 

f^d = ^ takes its maximal value [this means the inflection point of the function 



B = B{H)]. For H < a, the moment (dB)'^ is positive while for if > cr it is negative. 



2. The fourth order moment (SB)^ give informations about the deviation from the Gaus- 
sian approximation. This because it is implied in the so called excess coefficient (p3|). 
In the here discussed case we have: 

^ _kT M,H (S-3)e^p(-g) 
^ V /ioa4|A^exp(-f|) + l}' 



In probabilistic terminology |T3[ the distribution of a random variable is called lep- 
tokurtic, mesokurtic or platykurtic as the excess coefficient Ce satisfy the conditions 
Ce > 0, Ce = and Ce < respectively. Then by taking into account the expression 
of C^; one can say that for the here studied situation the fluctuations of the mag- 
netic induction B are leptokurtic, mesokurtic and platykurtic as the magnetic fleld 
strength H satisfy the conditions H > VSa, H = ^/Sa and H < ^/Sa respectively. 



IV. SUMMARY AND CONCLUSIONS 



1. We investigated the higher order moments of the fluctuations for complex thermody- 
namic systems (i.e. systems considered in the presence of flelds). Our approach uses 
a generalized statistical ensemble. We considered the case when the energy E, the 
usual extensive parameters {Xi}^^-^ and the fleld parameters {Yj}J^^ are fluctuating 
random variables. We flnd that the higher order moments of fluctuations depend on 
the fleld constraints. 
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2. The general results from the Sec. II were particularized for the case of a system situated 
in a magnetoquasistatic field. If the magnetic characteristics of such system are linear 
the third order moment of the magnetic induction and the excess coefficient are null. 
Therefore the description of the fluctuations for the magnetic induction can be done 
in the framework of Gaussian approximation. 

3. For nonlinear magnetic media (55)^ 7^ and consequently the fluctuations of B devi- 
ate from the normal distribution. Additionally the respective deviation are character- 
ized by the various value of the excess coefficient Ce given by the formula (^31). From 
the respective formula it results that the fluctuations of B can be leptokurtic (for 
H > VSor), mesokurtic (for H = VSa) and respectively platykurtic (for H < VScr) 
for the cases in which the function B = B[H) is given by the relation (|39D. 

Acknowledgments 

Our researches reported both here and in a previous paper were stimulated by the publi- 
cations of Professor Y. Zimmels. We express our gratitude to Prof. Zimmels for putting at 
our disposal a lot of his publications. We mention also that the present work takes advantage 
of some facilities of a grant from the Romanian Ministry of Education and Research. 



[1] S. Dumitru, A. Boer, Phys. Rev. E 64, 021108 (2001). 

[2] G. A. Korn, T. M. Korn, Mathematical Handbook Mc Graw Hill, New York 1968 (Russian 

version: Nauka, Moscow 1977). 
[3] S. Dumitru, Phys. Lett. 35A, 78 (1971). 
[4] S. Dumitru, Phys. Lett. 41A, 321 (1972). 
[5] S. Dumitru, Acta Physica Polonica A46, 149 (1974). 
[6] S. Dumitru, Optik 110, 110 (1999). 

[7] A. Miinster, Statistical Thermodynamics^ Springer- Verlag 1969. 

[8] Yu. B. Rumer, M. Sh. Ryvkin, Thermodynamics, Statistical Physics and Kinetics, Mir Pub- 
lishers, Moscow 1980. 
[9] G. Ruppeiner, Rev. Mod. Phys. 67, 605 (1995). 



10 



[10] Y. Zimmels, Phys. Rev. E 52, 1452 (1995). 

[11] Y. Zimmels, Phys. Rev. E 53, 3173 (1996). 

[12] Y. Zimmels, J. Chem. Soc, Faraday Trans., 94, 541 (1998). 

[13] [http://mathworld.wolfram.com/Kurtosis.html 



11 



